Double exchange magnets: Spin-dynamics in the paramagnetic phase 
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The electronic structure of perovskite manganese oxides is investigated in terms of a Kondo lattice 
model with ferromagnetic Hund coupling and antiferromagnetic exchange between f2 ff -spins using a 
finite temperature diagonalization technique. Results for the dynamic structure factor are consistent 
with recent neutron scattering experiments for the bilayer manganite Lai.2Sri.gMn207 . The sus- 
ceptibility shows Curie- Weiss behaviour and is used to derive a phase diagram. In the paramagnetic 
phase carriers are characterized as ferromagnetic polarons in an antiferromagnetic spin liquid. 
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Double exchange (DE) as an explanation of ferromag- 
netic order in doped manganites |l ] was proposed shortly 
after the first experimental work on these compounds 
H. Open questions are connected with the proper un- 
derstanding of the paramagnetic (PM) phase above the 
Curie temperature Tc and the giant magnetoresistance 
at the phase transition. Millis et al || argued that DE 
is neither sufficient to explain the size of the resistiv- 
ity above Tc nor the activated form of p(T) observed 
in many compounds, and suggested that Jahn- Teller po- 
laron formation is essential. A different point of view 
PL is that these systems are governed by strong corre- 
lations while spin disorder scattering is crucial for the 
understanding of the resistivity. Further complications 
not considered here may arise from orbital degeneracy 
H[| and charge disorder 0. 

A detailed investigation of the spin-dynamics is there- 
fore important for a deeper understanding of the mangan- 
ites. Inelastic neutron scattering experiments in the PM- 
phase were performed recently by Perring et al. Q for the 
highly doped bilayer manganite system Lai^Sri.gM^Oy 
H . Particularly surprising was the observation of grow- 
ing antiferromagnetic (AF) scattering at low energy par- 
allel to increasing ferromagnetic (FM) critical scattering 
as the FM Tc was approached from above. Our aim here 
is to shed light on these experiments by an investigation 
of the spin-dynamics within a standard model for these 
compounds. 

The appearance of ferromagnetism was explained by 
Zener jjj in terms of a strong Hund coupling Jh of the 
e g conduction electrons to S=3/2 core spins S? formed 
by the energetically well separated i2g-electrons. The e g 
electron in Mn 3+ couples to a S=2 high-spin configura- 
tion. A relevant model is the ferromagnetic Kondo lattice 
model ITcHjI 
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spins, respectively, while n- la = d\ a di a denotes the e g - 
density operator. A typical value for the Hund cou- 
pling is Jh ~ leV, while the nearest neighbor hopping 
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The operators d- ia and S? refer to e 9 -electrons and ti g 



0.15eV fgjbj. In view of the large local repul- 
' 8 eV, we shall furtheron exclude double occu- 
pancy of e g -orbitals, and define Jn/t = 6. The model 
explains ferromagnetism via the optimization of the ki- 
netic energy of the N e = n\a conduction electrons in 
the partially filled e g band when the spins are aligned 
at low temperature. At low doping manganite materi- 
als are usually AF insulators. Antiferromagnetic super- 
exchange interaction Jaf between the ^-electrons or 
spins is an obvious source of such correlations. 

The interplay of the global AF exchange interaction 
with the local FM alignment induced by the motion of 
e g -electrons via the DE mechanism is a subtle problem. 
A simple physical picture, which is consistent with our 
numerical study, is that of ferromagnetic spin polarons, 
where a local FM surrounding of t2g-spins moves with 
the holes through an otherwise AF background. 

Our calculations are based on a finite tempera- 
ture diagonalization approach developed by Jaklic and 
Prelovsek JT3] , which was successfully used for the study 
of the t-J model. Here due to the much larger Hilbert 
space only relatively small systems can be solved. Never- 
theless due to the high density of low energy excitations 
we believe that our results reflect already the physics of 
large systems. We have found large similarity between 
ID and 2D systems, in striking contrast to the physics of 
the t-J model. We shall present here mainly ID results, 
since the resolution in momentum space is larger. Our 
energy and temperature units are t = 1, while S = 1 is 
used instead of S = 3/2 i 2 g-spins to save computer time. 

A global measure for the evolution of magnetic corre- 
lations is given by the square of the total spin S = ^2 Sij 
which includes both core- and e 9 -spins. Results for dif- 
ferent e s -concentration and Jaf are shown in Fig. 1(a). 
The high-T limit is determined by independent spins (i.e. 
for T < Jh; for even higher T also the e g -spins decou- 
ple). In the absence of Jaf FM correlations develop as 
the temperature is lowered. For quite small Jaf this be- 
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haviour is already changed. For Jaf = 0.1 (S 2 ) decreases 
with decreasing temperature for a single hole. This sug- 
gests the evolution of AF correlations, while at higher 
doping FM correlations are established via the DE mech- 
anism. Further insight is gained from the uniform sus- 
ceptiblity X = (S^)/Nk B T. The corresponding l/x(T)- 
data, Fig. 1(b), show Curie- Weiss behavior, similar to 
experiments |l5|,[l(| . It is quite remarkable that for these 
small systems with N = 6 and 8 sites 1/%(T) follows a 
linear Curie- Weiss law to quite low temperatures, which 
allows us to determine a mean-field transition temper- 
ature Tg IF . Deviations at lower temperatures are ex- 
pected from non mean-field behavior of low-dimensional 
systems and finite size effects. 
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FIG. 1. Square of total spin (a) and inverse sus- 
ceptibility (b) vs. temperature for hole concentrations 
x = (l-N e JN) = 1/6, 3/6 (ID) and 1/8 (2D). A comparison 
with the FM polaron approach is shown for a single hole in 
ID (triangles) and 2D (squares) for Jaf ~ 0. 

The mean-field phase diagram (Fig. 2) reflects the in- 
terplay between antiferromagnetism and carrier induced 
ferromagnetism at higher doping. T^ IF scales with coor- 
dination number and is consistent with the experimental 
Tc-scale. The value chosen for Jaf{~ 0.1) is typical for 
real systems [||,|| . For Jaf = our data is well described 
by T£ IF ~ x(l - x) §J, while it differs in that respect 
from high-temperature series expansion results |]l7f . 

A simple estimate of (S 2 ) in terms of a ferromagnetic 
spin-polaron picture has been given by Varma (J] for the 
model without AF-interactions. The spin polaron is as- 
sumed to be determined by a number P of ferromag- 



netically aligned spins around the hole which allows the 
particle to improve its kinetic energy. The gain in kinetic 
energy is counterbalanced by the loss of spin entropy due 
to the coupling of P spins, which can no longer rotate 
independently. An estimate for the change of free energy 
is SE P = at/P r i + Pk B T\n(2S 2 + 1). The kinetic energy 
exponent r/ = 2/3, 1, 2 in 3, 2, and 1-dimcnsions, respec- 
tively. Minimization leads to P{T) = (^y) 1 ^ 1 ^ with 
B = 111(252 + !)■ Under the same assumptions the square 
of the total spin is given in terms of Sp = (S± + PS2) 



(S 2 ) = N 



cSp(S P + 1) + (1 - x - Px)S 2 (S 2 + 1) , (2) 



where in the case of few e g -holes x is the hole concentra- 
tion. Here S\ — 1 and S 2 — 3/2 (instead of 3/2 and 2). 
A comparison of this approach with exact diagonaliza- 
tion is shown in Figs. 1(a) and (b) for a ID (2D) system 
with 6(8)-sites and one e g hole. Good fits over the whole 
T-range were obtained by choosing smaller values for the 
exponent 77 = 1.5(0.6) for ID (2D), respectively. 

The effective spin S e f / for a single polaron may be de- 
fined as S e// (S e// + 1) = (S 2 ) T - (S 2 ^. Near Tg F 
we obtain S e ff ~ 7 for the two-dimensional case. This 
corresponds to an enhancement factor ~ 14 with respect 
to the spin S = 1/2 of the e g conduction electron (hole). 
This result is qualitatively consistent with EPR exper- 
iments [M where large effective spins at temperatures 
slightly above Tc were observed. 
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FIG. 2. Mean-field phase diagram of a ID chain for 
Jaf = 0.1 (triangles) vs. (1 - x) = N e /N. The bound- 
aries between AF, FM and PM phases (solid lines) are guides 
to the eye. (Anti)periodic boundary conditions were used for 
odd (even) electron number [18]. The Jaf = data (squares) 
is compared with the predicted behavior Tq IF ~ x(l — x). 

The study of correlation functions (CF) gives further 
support to the FM polaron picture even at high dop- 
ing concentration and with frustrating AF interactions, 
although the simple formulation above is certainly not 
sufficient for these cases. The temperature dependence 
of the (Si • Sj) correlation functions for Jaf — 0.1 is 
shown in Fig. 3. We distinguish (a) a low- and (b) a 
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high-doping case where the system approaches the AF 
(FM) phase (Fig. 2) at low temperature, respectively. In 
the single hole case (a) the nearest-neighbor (n.n.) CF 
becomes AF as T is lowered, while further neighbor cor- 
relations remain small. The CF's for x = 0.5 (b) become 
all FM at low T. Yet it is interesting to notice that for 
intermediate T the n.n. CF is antiferromagnetic. For 
T > t all CF's vanish. 

To investigate the structure of the spin-polaron we cal- 
culate the relative correlation function (RCF) 



C" +1 = 



(n^Sj+n • Si+n+l), 



(3) 



which measures the spin-CF (So • Si) in a distance n from 
the position of the moving hole. Here n\ = (1 — n{) and 
1 is a lattice vector. The results in Fig. 3 show the FM 
alignement of the spins at the position of the hole and on 
the neighboring site, i.e. Cq, at low T even in the single 
hole case where (So -Si) becomes antiferromagnetic. The 
large drop at very small T suggests that FM-polarons 
are unstable in the low temperature AF-phase consistent 
with Fig. 1(a). The CF's Cq are weakly antiferromag- 
netic at higher temperatures similar to (Si • Sj). 




FIG. 3. Temperature dependence of spin correlation func- 
tions (SiSj) and C„ +1 (Eq. 3) for one (a) and 3 e s -holes (b) 
for a 6-site chain with Jaf =0.1. 



In the single hole case (a) all other (n ^ 0) RCF's 
become gradually AF with decreasing temperature, i.e. 
indicating that the AF-phase dominates. The CF's for 
x = 0.5 (b) are indicative for frustrated antiferromag- 
netism at temperatures above T ~ 0.5, while at lower 
temperature all CF's become ferromagnetic. Hence on 
the basis of these equal-time CF's there is no evidence for 
significant AF correlations at low T near T(f F for high 



doping. Nevertheless there are clear AF-correlations in 
the low-frequency response. 

The complete information about the spin-dynamics is 
contained in the dynamic structure factor 
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Sa(u) = — dtexp (-iujt)(SqS-q(t)), 

27r J-oo 



(4) 



which can be directly measured by inelastic neutron scat- 
tering. A central result of this work is the evolution of 
Sq(ui) with temperature (Fig. 4). As an example results 
for a 6-site chain with 3 e g -holes are shown for Jaf = 
(a) and 0.1 (b) at different momenta. At high T critical 
scattering extends over the whole Brillouin zone, while at 
low temperature it is confined to q ~ in the absence of 
AF-interactions. In case (a) FM magnons emerge when 
lowering the temperature already well above T^ F ~ 0.1. 
Their energy scale is consistent with theoretical work on 
the magnon dispersion |l9[ ] and neutron scattering in the 
FM-phase |2(J of La .7Pb .3MnO3. Antiferromagnetic 
frustration appears to be small in this system, since al- 
ready LaMnOa consists of FM planes with a well devel- 
oped magnon spectrum even at Tc [ plf . The interlayer 
coupling is AF and weak in this system. 
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FIG. 4. IVlagnetic fluctuation spectrum S(q, u) as function 
of oj for different momenta and temperatures. Data for a 6-site 
chain with e 9 -electron concentration x = 0.5 for (a) Jaf = 
and (b) Jaf = 0.1. 

In the frustrated case (b) the magnon structures are 
hardly visible, consistent with a substantial softening. 
The dominant feature is critical scattering for all q. As 
temperature is lowered FM-critical scattering is dom- 
inant, yet surprisingly low energy AF-correlations at 
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q = 7r are growing as well. We note that our data is 
broadened by 6 = 0.05. The natural linewidth of the 
lo r*> structures appears to be smaller than T. 

The low frequency response <5<S'(q) = J _° dojSq(u}) in- 
tegrated over a small energy window (u> c = 0.05), is fre- 
quently measured in neutron scattering experiments 
It is interesting to contrast the momentum dependence 
of SS(q) from the behavior of the static structure factor 
5(q) = /_ d,LuSq(u>) , which measures equal-time corre- 
lations (Fig. 5). In the high-temperature limit 5(q) is 
constant, reflecting the absence of spin-correlations. At 
low temperature S (q) develops a peak centered at q = 
in the model without antiferromagnetic interactions due 
to carrier induced ferromagnetic correlations. With in- 
creasing Jaf the data shows the evolution of short-range 
(nearest-neighbor) AF correlations. 
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FIG. 5. Comparison of S'(q) and integrated low- frequency 
response SS(q); parameters as in Fig. 4. 

The observation of pronounced AF-correlations in 
<55(q), e.g. for Jaf — 0.1, in the presence of FM criti- 
cal scattering is consistent with the experimental finding 
of local AF correlations by Perring et al ||. These ex- 
periments show strong AF correlations in the PM-phasc 
(at low frequency) and their sudden disappearence at Tq 
when the low-T FM-order is established. We stress that 
the observation of a well defined AF peak in the low- 
frequency spin- fluctuation spectrum Sq(u> — > 0) is ac- 
tually very surprising since these experiments were per- 
formed for a highly doped (x — 0.4) double-layer system 
and at temperatures slightly above the transition into 
the FM-phase. We find considerable qualitative agree- 
ment between the present model and neutron scattering, 



although we do not find a significant peak structure in 
the momentum dependence at q — ir, except for larger 
values for Jaf where the FM scattering is reduced. The 
disappearance of the AF signal in the low-temperature 
FM phase |8) is consistent with our studies p^ |. 

In summary we have shown, that the competition 
between AF and FM correlations detected by neutron 
scattering in the highly doped double layer compound 
Lai.2Sri.gMn207 is consistent with the physics of the 
frustrated Kondo lattice model. From our study of vari- 
ous correlation functions, we conclude that the tempera- 
ture and doping dependence of the magnetic correlations 
in the paramagnetic phase can be qualitatively under- 
stood in terms of FM polarons moving in a short-range 
AF-spin liquid. Slow spin-dynamics as a consequence of 
frustration may lead to localization of magnetic polarons. 
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